Non-adiabatic control of quantum energy transport in ordered arrays 
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An elementary excitation in an aggregate of coupled particles generates a collective excited state. 
We show that the dynamics of these excitations can be controlled by applying an instantaneous 
external potential that modifies the quantum phase of individual particles. We show that spe- 
cific phase transformations can be used to accelerate or decelerate quantum energy transport and 
spatially focus delocalized excitations onto different parts of ordered arrays of quantum particles. 

PACS numbers: 42. 50. Hz, 34.50. Ez, 42.50.Dv, 37.10. Pq 



Many physical systems, such as molecular dye 
aggregates pQ, photosynthetic complexes [2], atoms and 
molecules in optical lattice traps [3] H] or self-assembled 
monolayers on solid substrates [5], represent a chain of 
coupled monomers, which, once excited, can exchange 
the excitation energy. This leads to energy transport. 
Understanding and controlling energy transport is cru- 
cial for many applications such as building efficient light- 
harvesting devices [6] and scalable quantum computing 
networks [7]. In this work, we explore the general algo- 
rithm for controlling transport through an array of cou- 
pled quantum monomers by applying monomer-specific 
external perturbations. This is possible to realize in an 
ensemble of spatially separated monomers, such as Ryd- 
berg atoms or molecules trapped on an optical lattice. 

Excitation of a coupled many-body system generates 
a wavepacket representing a coherent superposition of 
single-site excitations. We explore the possibility of 
using external fields to control their dynamics. The 
method proposed here is based on shaping the many- 
body wavepackets by a series of sudden perturbations, 
in analogy with the techniques developed for strong-field 
alignment and orientation of molecules in the gas phase 
[8] . Alignment is used in molecular imaging experiments 
and molecular optics [8[ [9] , and is predicted to provide 
control over mechanical properties of molecular scatter- 
ing [10 . Our technique allows for the possibility to study 
condensed-matter excitations and energy transport with- 
out statistical averaging. We demonstrate that a de- 
localized excitation can be spatially focussed, enabling 
directed transport of energy to a particular monomer. 
This can be exploited for the experimental realization of 
controlled chemical transformations and as a toolbox for 
optimizing photon-harvesting systems. 

We consider an ensemble of N coupled identical 
monomers possessing two internal states \g) and |e) ar- 
ranged in a one-dimensional array with translational 
symmetry. The excitation spectrum of the system is de- 
termined by two parameters: the monomer excitation 
energy AE e _ g and the coupling strength a(n) between 
the monomers separated by the distance na, where a is 
the nearest neighbor separation and n is an integer. The 



singly excited state of the system is 

N 

|^xc) = Ec„|e„)nift), (1) 

n=l iy^n 

where \g n ) and \e n ) denote the ground and excited states 
in site n, C n = e iakn /y/N and |^ exc ) rep- 
resents a quasi-particle called Frenkel exciton, charac- 
terized by the wave vector k [11 j. The energy of the 
exciton is given by E(k) = AE e - g + a(k) with a(k) = 
^2 n a(n)e~' iakn . In the nearest neighbor approximation, 

E(k) = AE e _ g + 2a cos ak, (2) 

where a = a(l). A localized excitation placed on a single 
site (or a limited number of M < N sites) is a coherent 
superposition of the exciton states IVwW) with differ- 
ent k. 

Control over energy transport in an ordered array can 
be achieved by (i) shifting the exciton wavepackets in the 
momentum representation (which modifies the group ve- 
locity and the shape evolution of the wavepackets) and 
(ii) focusing the wavepackets in the coordinate represen- 
tation to produce localized excitations in an arbitrary 
part of the lattice. To achieve this, we propose to apply 
a series of site-dependent perturbations. For (i), the es- 
sential idea is to add a factor e %5an to each term in the 
expansion ([l]), so that each |Vw(&)) is transformed into 
|V ; exc(& + ^)/^ This transformation shifts the wavepackets 
by 5 in /c-space while preserving their shape. As a result, 
one can engineer wavepackets probing any part of the 
dispersion E(k) leading to different group velocity and 
shape evolution. 

Adding a site-dependent phase to the excitonic wave- 
function exploits an interplay of the adiabatic and sud- 
den time scales. Consider, first, the n-th monomer sub- 
jected to an external field £ n (t) which varies from to 
some value and then back to in time T. If the varia- 
tion is adiabatic with respect to the evolution of the free 
monomer states, T ^> h/AE e _ g: each eigenstate |/) of 
the monomer acquires a state-dependent phase shift [12] 

|/„(T)> = e-^|/„(0)>, where <p{ = ± [ E*(t)dt, (3) 

(t) is the instantaneous eigenenergy and / can be e or 
g. Now consider the action of such phase change on the 
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collective excitation state ([I]). If T <C h/a, the change 
is sudden with respect to the excitation transfer between 
monomers and the state ([I]) acquires a site-dependent 
phase $ n = (\f n — (\> 9 n . If $ n = $o + na5, then the mo- 
mentum 5 is imparted onto the excitonic wavefunction. 
By analogy with "pulsed alignment of molecules" [8] , we 
call this transformation a "phase kick" or "momentum 
kick" . Its action is also similar to that of a thin prism on 
a wavefront of a monochromatic laser beam. 

In order to illustrate the shifting of exciton wavepack- 
ets in the momentum space, we solve numerically 
the time-dependent Schrodinger equation with the un- 
perturbed Hamiltonian H exc = AE e _ g ^ n \e n )(e n \ + 
^Cn,m a ( n ~ m )\ e n,gm)(gn,e m \, subjected to a transient 
site-dependent external perturbation. The excitation at 
t = is described by a Gaussian wave packet of the exci- 
ton states \i/j exc (k)} : with the central wavevector k = 0. 
The external perturbation is provided by a spatially in- 
homogeneous pulse of electromagnetic field. Fig. [I] shows 
that the entire wave packet acquires momentum during 
the pulse (left panels). This is manifested as a phase 
variation in the coordinate representation, and as a shift 
of the central momentum in the /c-representation. After 
the external perturbation is gone, the wave packet does 
not evolve in the ^-representation and moves with the ac- 
quired uniform velocity in the coordinate representation. 

In optics, a thin lens focuses a collimated light beam 
by shifting the phase of the wavefront, thus converting 
a plane wave to a converging spherical wave. Similarly, 
a phase kick can serve as a time domain "lens" for col- 
lective excitations: an excitation delocalized over a large 
number of monomers can be focused onto a small re- 
gion of the array. By analogy with optics, a concave 
symmetric site-dependent phase <£(n) applied simultane- 
ously to all monomers may turn a broad initial distribu- 
tion C n (t = 0) into a narrow one. The dynamics of the 
amplitudes in Eq. ([I]) is given by 

C m (t) = ^^2C n (t = Oymn)+ka(m-n)-E(k)t/h] ^ 
n,k 

The shape of the dispersion curve E(k) is crucial for the 
focusing efficiency: while waves with quadratic dispersion 
can be focused, those with linear dispersion propagate 
without changing the wavepacket shape [13] . The cosine- 
like dispersion of excitons ([2| offers both the quadratic 
(at low k) and linear (at k ~ ±7r/2a) regions. A rigorous 
analysis shows that, for the cosine dispersion, the wave 
packet (l) with C n = J n _ no (2ar)e i7r ( n - n °)/ 2 , where J n is 
the Bessel function of order n, must focus onto the single 
site no at time r. However, creating such a wavepacket 
may require multiple and complicated phase transforma- 
tions. For practical applications, it may be easier to con- 
strain all phase transformations such that the wavepacket 
always remains in the quadratic part of the cosine dis- 
persion pi. For such wavepackets, adding a quadratic 
phase <£(n) = &o(n — no) 2 leads to focusing around site 
no- Below we illustrate the effect of the quadratic phase 
transformation for two types of initial states. 
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FIG. 1: (Color online) Evolution of the exciton wavepacket in 
the k and real spaces. The phase of the wave function is shown 
by color. The calculation is for a one-dimensional array of 
201 monomers with a = 22.83 kHz and AE e - g = 12.14 GHz, 
which corresponds to LiCs molecules trapped on an optical 
lattice with a = 400 nm and subjected to a homogenous DC 
field of 1 kV/cm directed perpendicular to the intermolecular 
axis. The kicking potential is provided by a A = 1024 nm 
Gaussian laser beam, with the propagation direction along 
the array axis, focused to 5 /im, with the intensity at the 
focus equal to 10 7 W/cm 2 . The laser pulse is on between 
and 3 /is. The molecules are placed on the beam axis with 
the first molecule 5 fim away from the focus. Thus, the laser 
intensity is gradually decreased along the array. 



First, consider a broad Gaussian wavepacket ([I]) with 
C n (a x ;t = 0) = \/a/cr x ^ exp [-a 2 (n - n ) 2 /2a x ] ] 
where <j x ^> a. By applying the quadratic phase to the 
Gaussian wave packet with quadratic dispersion, one cre- 
ates a wave packet with the width in the /^-representation 

m 

tr k {<r x , *o) = — v/l + 4$§<74/ ffl 4. (5) 

Eq. ([5| shows that the value of $0 determines the de- 
gree of focusing. However, for systems with dispersion 
([2]), large values of $0 may take the wavepacket out- 
side the quadratic part of the dispersion. To find the 
optimal phase $q ^ na ^ keeps the wavepacket within the 
quadratic dispersion while focusing it, we use the con- 
dition Afc = acjk ^ 1> which yields $0 = a/2<j x for the 
optimal focusing. At time t* ~ l/4a$Q, the wavepacket 
is most focused and has a width 

a x F (*S) = ° x = w a. (6) 
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Second, consider a completely delocalized excitation 
jlj with C n (k;t = 0) = e iakn /VN describing an eigen- 
state of an ideal system of N coupled monomers. If E(k) 
in Eq. Q is assumed to be E(k) = AE e - g — aa 2 k 2 , 
the phase $(n) = § (n - n ) 2 with $ = tt/N re- 
duces the amplitudes C m at t* = N/4tt to C m (k;t*) « 
Vie- iiVa2fc2 / 47r ^ n , no _^, where i/ fe = Nak/2ir [14]. How- 
ever, an estimate shows [14] that this phase modulation 
would create a wavepacket with the width Ak(®o) ~ 
2N&o in the momentum space spanning the entire range 
of k- values. For an arbitrary value of Afc(<I>o), the site 
amplitudes at the time of focusing t* = l/4a<I>o are [T4] 



p iA($ )n 2 /2^ 



2z 



C n (k-t = U) * J /t n sin(nA fc ($ )/2). 

n y 7rA fe ($ ) 

(7) 

In order to keep the linear part of the dispersion spectrum 
unpopulated, we choose $q ~ l/2iV so that A fc ($o) ~ 1- 
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FIG. 2: Focusing of a completely delocalized collective exci- 
tation (panels a and b) and a broad Gaussian wave packet of 
Frenkel excitons (panels c and d) using the quadratic phase 
transformations at t = as described in text. The dashed 
lines show the initial distribution magnified by 20 and 5 re- 
spectively in (b) and (d). The solid curves in panels (b) and 
(d) correspond to two different phase transformation focus- 
ing the same wave packet onto different parts of the array. 
The calculations are performed with the same parameters a, 
a, and AE e - g as in Figure 1. The results in panels (a) and 
(b) are computed with nearest and second neighbor couplings 
only; the results in panels (c) and (d) are computed with all 
couplings accounted for. 

Eqs. (6) and (7) are valid for a many-body system 
with nearest neighbor interactions only. In most physical 
systems, the energy dispersion is modified by long-range 
couplings. In order to confirm that the above predictions 
are also valid for systems with long-range interactions 
and illustrate the focusing of delocalized excitations, we 
compute the time evolution of the wave packets by solv- 



ing the wave equation numerically for a system with long- 
range dipole - dipole interactions. Figure 2 illustrates the 
focusing dynamics of a completely delocalized excitation 
in a system with nearest neighbor and second-order cou- 
plings (panels a and b) and a broad Gaussian wave packet 
in a system with all long-range couplings explicitly in- 
cluded in the calculation (panels c and d). The results 
show that, although the long-range couplings lead to an 
oscillating interference pattern, the collective excitations 
can be focussed to a few lattice sites. 

By analogy with optics, it should be expected that the 
focusing phase transformations must be more effective in 
systems of higher dimensionality. To illustrate this, we 
repeated the calculations presented in Figures 2c and 2d 
for a delocalized excitation placed in a 2D array with an 
external potential that modulates the phase as a func- 
tion of both x and y. Figure 3 shows the focusing of an 
initially broad wave packet onto different parts of a 2D 
lattice induced by the quadratic phase transformation. 
The calculations include all long-range couplings. Fur- 
ther, similar to spatial light modulators in optics, multi- 
ple phase kicks of varying profiles separated by periods 
of free evolution may enable arbitrary shaping of the ex- 
citonic wave packet in the coordinate and ^-spaces. 
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FIG. 3: Focusing of a delocalized excitation in a 2D array 
shown at t = in panel (a) onto different parts of the lat- 
tice (panels b-d). The probability distribution in panel (a) 
is enhanced by the factor 72 for clarity. The calculations are 
performed with the same parameters a, a, and AE e - g as in 
Figure 1 and the quadratic phase transformation at t — 0. 



While the results shown in Figures 1-3 are general, for 
concreteness, the calculations presented are performed 
for an experimentally realizable system of polar diatomic 
molecules in the 1 E electronic state trapped on an op- 
tical lattice [15]. We assume that the optical lattice 
is prepared with one molecule per lattice site and that 
tunneling of molecules between different lattice sites is 
suppressed. The rotational states of the molecules are 
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labeled by the quantum number of the rotational angu- 
lar momentum J and the projection Mj of J on the 
space-fixed quantization axis Z. We choose the rota- 
tional ground state \J — 0,Mj = 0) as \g) and the 
rotational excited state \J = l,Mj = 0) as |e). The 
state | J = 1, Mj = 0) is degenerate with the states | J = 
l,Mj = ±1). This degeneracy can be lifted by apply- 
ing a homogeneous DC electric field, making the \g) and 
\e) states an isolated two- level system. The molecules in 
different lattice sites are coupled by the dipole - dipole 
interaction V&&(n — m). The separation between adjacent 
sites of an optical lattice is equal to half the wavelength 
of the trapping laser and is, typically, in the range be- 
tween 256 and 800 nm [4 . The magnitude of the coupling 
constant a(n — rn) = (e ni g m \V&&(n — m)\g ni e m ) between 
molecules with the dipole moment 1 Debye separated by 
500 nm is on the order of 1 kHz [TB] . 

For molecules on an optical lattice, one can implement 
the phase kicks by modifying the molecular energy levels 
with pulsed AC or DC electric fields. The rotational 
energy levels for 1 £ molecules in a combination of weak 
AC and DC electric fields are given by [T7] 



E 



J,Mj 



BJ{J+l) 
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^G(J,Mj) 



a ± £ 2 



AC 



( a \\ 



oli_)£\c 



F(J,Mj) (8) 



where B is the rotational constant, G(0, 0) = —1/3, 
G(1,0) = 1/5, F(0,0) = 1/3, G(1,0) = -1/5, £ AC is 
the envelope of the quickly oscillating AC field, an and 
a± are the parallel and perpendicular polarizabilities and 
H is the permanent dipole moment of the molecule. 

The momentum shift of the exciton wave packets can 
be achieved by applying a time-varying DC electric field 
£(t) = £* + £ (n) sin 2 (7rt/T), where £{n) is linear with 
respect to n. Assuming that £ (n) = (n — no) A and 
£ (n) <C £*, and using Eqs. ^ and ([8]), gives S = 
4A£*[i 2 T/15hBa. We have confirmed this result by a 
numerical computation showing that for LiCs molecules 
in an electric field of £* = 1 kV/cm, an electric field pulse 
with A = 7.434 x 10 -4 kV/cm and T — 1 /is results in a 
kick of 5 = 7r/2a, bringing an excitonic wave packet from 
the k = region to the middle of the dispersion zone. 

An alternative strategy is to use a pulse of an off- 
resonance laser field. Consider a laser pulse whose in- 
tensity varies linearly across the molecular array: 



[£ 2 + (n - no)^ 2 ] sm 2 (irt/T) (0 < t < T) . 

(9) 

Using Eqs. ^ and (J8|, we estimate the momentum kick 
by such a pulse as S = T£ 2 (a\\ — a±) /15a. The results 
shown in Figure 1 correspond to this set of parameters. 

In addition to the momentum kicking technique, en- 
ergy transport in an array of molecules on an optical lat- 
tice can be controlled also by varying the angle between 
the intermolecular axis and the applied DC field. This 
angle influences the dipole-dipole interaction between the 
molecules and determines the sign and magnitude of a, 
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and therefore the shape of the dispersion curve, as shown 
in Figure 4(a). This enables control over the sign and 
magnitude of the group velocity of an excitonic wave 
packet with k ^ 0. Dynamically tuning one can propa- 
gate a localized excitation to different parts of the lattice, 
as shown in Figure 4b. 
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FIG. 4: (Color online) (a) Exciton dispersion curves for a 
ID ensemble of diatomic molecules on an optical lattice for 
different angles between the direction of the external DC 
electric field and the axis of the molecular array; (b) Propa- 
gation of a wave packet centered at ak — —1.0464 controlled 
by tuning the electric field direction. Thin dashed line depicts 
the corresponding angle variations with time. 

In summary, any elementary excitation in an aggregate 
of coupled monomers can be represented as a coherent su- 
perposition of free exciton states. We propose a method 
for shaping the exciton wave packets using "phase kick- 
ing" perturbations that temporarily modulate the energy 
levels of the monomers leading to a displacement of the 
wave packets in the momentum representation. This can 
be used to control the energy transport in ordered en- 
sembles of quantum particles. In particular, we show 
that the phase kicking technique can be used for focus- 
ing a delocalized collective excitation to a particular part 
of a ID or 2D lattice. 

Control over excitation transport is needed for creat- 
ing networks of quantum processors where information is 
transmitted over large distances with photons and stored 
in arrays of quantum monomers via one of the quantum 
memory protocols [18]. Momentum kicking can be used 
for information transport within a single array. Focusing 
excitonic wave packets enables local storage of informa- 
tion, while controlled interactions of multiple excitons 
[T9] or excitons with lattice impurities [20 may be used 
to implement logic gates. 

Controlled energy transport in molecular arrays may 
also be used for the study of controlled chemical inter- 
actions for a class of reactions stimulated by energy ex- 
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citation of the reactants. Directing quantum energy to Studies, 
a particular lattice site containing two or more reagents 
can be used to induce a chemical interaction [21], an in- 
elastic collision or predissociation [22] with the complete 
temporal and spatial control over the reaction process. 
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